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Low-temperature magnetization curves of the classical frustrated ferromagnetic chain in the ex- 
ternal magnetic field near the transition point between the ferromagnetic and the helical phases is 
studied. It is shown that the calculation of the partition function in the scaling limit reduces to the 
solution of the Schrodinger equation of the special form for the quantum particle. It is proposed that 
' the magnetization of the classical model in the ferromagnetic part of the phase diagram including 

^N) , the transition point defines the universal scaling function which is valid for quantum model as well. 

Explicit analytical formulae for the magnetization are given in the limiting cases of low and high 
O ' magnetic fields. The infiuence of the easy-axis anisotropy on the magnetic properties of the model is 

I studied. It is shown that even small anisotropy essentially changes the behavior of the susceptibility 

in the vicinity of the transition point. 

■ I. INTRODUCTION 

B' .... 

I ' Strongly frustrated low-dimensional magnets have attracted much attention last years 1 1 . A very interesting class 
, of such compounds is edge-sharing chains where CuO^ plaquets are coupled by their edges l-Q- An important feature 
of the edge-sharing chains is that the nearest- neighbor (NN) interaction Ji between Cu spins is ferromagnetic while 
the next-nearest-neighbor (NNN) interaction J2 is antiferromagnetic. The competition between them leads to the 
' frustration. A minimal model describing the magnetic properties of these cuprates is so-called F-AF spin chain model 
the Hamiltonian of which has the form 

i7 = Ji Vs„.S„+i + J2Vs„.S„+2-/iV5^ (1) 



o 



- '—I 

X 



where S„ is the spin operator on ri-th site, h is the external magnetic field and the exchange integrals are Ji < and 
J2 > 0. 

This model is characterized by the frustration parameter a — J2/ | ^1 1 ■ The ground state phase diagram of the 
quantum s = i model has been intensively studied [8l-[l6| . The ground state of model ([1]) at ft, = is ferromagnetic for 
^ I ■ a < i . At a = i the quantum phase transition to the phase with incommensurate spin correlations of the helical type 
\ takes place. Remarkably, this transition occurs at the same frustration parameter ct = j both in the quantum and in 
t — ■ the classical model. However, the influence of the frustration on the low-temperature thermodynamics in the vicinity 
' of the transition point is less studied. This problem is of a special interest because recently studied edge-sharing 
compound Li2ZrCuOi is well described by the F-AF model with the frustration parameter close to a = ^ 

At present the low-temperature thermodynamics of the quantum s = i model at a ^ can be studied only either 
by numerical calculation of finite chains or by approximate methods. On the other hand, the classical version of 
model ([TJ can be studied exactly at T — > and the classical limit is a starting point for the study of quantum 
effects. Another reason to study the classical version of F-AF model comes from the following argument established 
for the quantum spin-s ferromagnetic chain, i.e. for model ([T]) at a = 0. It was conjectured in Ref.fis'l that the 
^ , low-temperature magnetization of this model is a function of the scaling variable gp = \Ji \ h/T^. According to this 
I scaling hypothesis the normalized magnetization m — (S^) / s {{S^) is the magnetization per site) of the quantum 
chain is expressed as 

m{T,h) = ^{gF) (2) 

This equation is valid in the scaling limit, which means that T — and ft — but (7_f is fixed. Then the dependence 
of the magnetization m on the spin magnitude s comes only via the scaliiiR variable gp. Generally, the calculation of 
the function 4){gF) is a very complicated problem. It was proposed in Ref.[l8| that this function can be obtained from 
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the solution of the classical ferromagnetic chain and such scaling function 4>{gF) was obtained explicitly in Ref . [Tsl, \\i 
In particular, the zero-field susceptibility x is 



dm 
" 'dh 



Actually the conjecture of the universality of the function 4'{gF) is based on the following observations [lal20]: the 
zero-field susceptibility of the s — ^ Heisenberg ferromagnetic chain at T — > coincides with that given by Eq.([3]); the 
magnetization to(T, h) obtained numerically from the thermodynamic Bethe-ansatz equations and plotted as a function 
of qf = h/ST"^ approaches to 4>{gF) at T — 0; the leading terms of the spin-wave expansion for magnetization coincide 
with those for 4>{gF)- In addition, as noted in Ref.fl^, the hypothesis of the universality originates in the universal 
behavior of the spin- wave excitations from the ferromagnetic ground state for both quantum and classical model. For 
this reason it is naturally to expect that such universality remains for all a corresponding to the ferromagnetic ground 
state, i.e. for < a < |. Moreover, the function (pigF) for < a < -j (but a not too close to |) will be the same 
as for a = but with gF replaced by (7f(1 — 4a). Really, the zero- field susceptibility x = 2s* (1 — Aa)/'iT'^ fits very 
well with numerical and analytical results (21j . However, vanishes ai a = ^ signalling the change of the critical 
exponent at the transiti on p oint. 

In our previous paper [23] we studied the zero-field susceptibility of the classical F-AF chain exactly at the transition 
point a = ^ and we have shown that x ~ T~*/^ in contrast with the low-temperature asymptotic x ^ T^^ for 
< a < i. As was shown in Ref. [22| the change of the critical exponent for x is a consequence of the modification 
of the energy of spin-wave excitations from e{k) '-^ fc^ for < a < | to e{k) ~ fc** at a = |. Therefore, the form of 
the universal magnetization curve and the scaling variable for a = i (if the universality is valid) differ from the case 
< < -J and require a special study. 

Another interesting problem related to the F-AF model is the influence of the anisotropy of exchange interactions 
of the easy-axis type on the low-temperature magnetic properties of this model. This problem is actual because it 
is known that in the real ed ge-s haring compounds the exchange interactions are anisotropic and the anisotropy can 
be of the easy-axis type (23l. |24| . Though this anisotropy is weak it can be important especially near the transition 
point. In particular, it essentially changes the behavior of the zero-field susceptibility [25| . 

In this paper we investigate the effect of weak anisotropy on the magnetic curves of the classical F-AF model at 
the transition point. In the low-temperature limit the easy-axis anisotropy of the NN and NNN interactions have the 
same effect (we will explain this fact below) and for simplicity we consider the anisotropy of the NN interaction only, 
i.e. we add to Hamiltonian ([l} the term 

-(A -1)^5^5^+1 (4) 

where A > 1. 

It is interesting to note that for the pure ferromagnetic case (a — 0) the similarity in the magnetic properties of 
quantum and classical models remains in the case of the easy-axis anisotropy. This resemblance is based on the close 
relation between the classical solitons and the quantum multimagnon bound complexes. In this paper we will elucidate 
the question to which extend the resemblance between the anisotropic quantum and classical models remains in the 
F-AF model. 

The paper is organized as follows. In Section II the continuum version of the model is introduced and the scaling 
parameters are determined. The calculation of the partition function is reduced to the solution of the Schrodinger 
equation of a special type. In Section III the behavior of the magnetization curve at the transition point is studied. 
The asymptotics of magnetization for low and high magnetic field are presented and relation to the quantum spin 
model is discussed. The numerical and analytical results for the magnetization curve in the helical phase are given in 
Section IV. In Section V the infiuence of the easy-axis anisotropy on the magnetic properties is studied. The summary 
of the obtained results is given in Section VI. 



II. PARTITION FUNCTION IN THE CONTINUUM LIMIT 

In Refs.[22l. [26j we studied the partition function and the spin correlation functions of the classical F-AF chain in 
the vicinity of the transition point a = -j at zero magnetic field. This study was based on the use of a continuum 
approximation and the interpretation of the partition function as a path integral for the quantum particle in a potential 
well. However, the extension of the model to non-zero magnetic field and/or non-zero anisotropy needs the essential 
modification of this approach. 
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In the vicinity of the transition point a = i it is convenient to rewrite Hamiltonian ([T]) with the anisotropic term 
(HI in the form 



H=l E(Sn+i - 2S„ + S„_i)2 -ha- 5^(S„+2 - S„)2 - (A - 1) ^ S^S^^, - E ^« 



(5) 



In Eq.([5]) we put |Ji| — 1 and omit unessential constant. 

In the classical approximation the spin operators S; are replaced by the classical vectors Si — sfti, where Hi are 
the unit vectors. In the low-temperature limit the thermal fluctuations are weak so that neighbor spins are directed 
almost parallel to each other. Therefore, at T — > we can use the continuum approximation replacing fii by the 
classical unit vector field n{x)^ so that 



i+l 



215*^ + S'i_i) 

{Si+2 — Si) 



dn[x) 



2s 



dx 



(6) 



where the lattice constant is chosen as unit length. 

Using Eqs.® Hamiltonian ([5]) goes over into the energy functional 



E [n{x)] 



dx 



an 
dx"^ 



s2(4a - 1) ( dn 



— shuz 



(7) 



One can easily check that in the continuum approximation the easy-axis anisotropy of NNN interactions A2 results 
in the term as^ (A2 — l)"-^- This term merely changes the coefficient in the third term of Eq.(I7]), so that the results 
obtained below cover the case of the NNN anisotropy as well. 

Energy functional ([7]) contains the second order derivative d^n/dx'^ in contrast with that for the ferromagnetic chain 
[isj . which contains the first order derivative dn/dx only. This fact demonstrates an essential difference between the 
cases a = and a = t. 



The partition function is a functional integral over all configurations of the vector field on a ring of length L 

E [n{x)] ' 



It is useful to scale the spatial variable as 



D [n{x)\ exp 



2^1/3, 



T 



g2/3 



Then, the partition function takes the dimensionless form 



Z = j i?[n(0]expi-^' 



d^ 



where 



d^fi 

de 



,s2/3 



J f dn\ 



-5nt 



gnz 



is the scaled system length and 
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(4a - 1)5^/3 



(A - 1)58/3 



hs^ 

j-4/3 



(8) 



(9) 



(10) 



(11) 



(12) 



are the parameters of the model scaled by temperature. 

As follows from Eq. ()10p the partition function and with it the low-temperature thermodynamics of the F-AF model 
near the transition point is governed by three scaling parameters 7, 6 and g. The definition of the scaling parameters 
implies that we consider the scaling limit when T — 0, a — >■ A — 1, /i — >• 0, but the values of the scaling 
parameters 7, S and g are finite. 



We express the unit vector field through two scalar fields 0(^) and n^(^) 

n{0 = {cos(f>y/l - n2, sin 0^/1 - n^^n^) 



(13) 



and the magnetic field is directed along the Z axis. In terms of the fields ^iO the partition function takes 

the form of the functional integral 



(14) 



where the energy density W{nz, (j)) has a rather cumbersome form: 



W 



1 — n; 



2nznj.nz 



2nznz(t> 



1 — n 
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1 — n: 



+ (1 - n 



2\12 



gnz 



(15) 



Here hz, riz and (f>, 4> are the first and the second-order derivatives of Uz and (j) with respect to ^. 

If we treat ^ as an imaginary time then partition function (1141) takes the form of a path integral of a quantum 
particle with the Euclidean Lagrangian W{nz, (f>)- Here we notice that W{nz, 4>) comprises the derivatives (j) and </> of 
the field (j>{£,), but does not contain explicitly the field 0(f) itself. This allows us to rewrite partition function ([Ti)) in 
terms of a new field 



(16) 



The energy density W{nz,4>) contains explicitly the field nz(f) and its derivatives hz and hz- Presence of the 
second-order derivative requires the use of the special methodology developed by Ostrogradski [131 which allows to 
obtain the Hamiltonian corresponding to the higher gradient Lagrangian. In the Ostrogradski formalism [2^ the 
independent generalized coordinates are and v = hz- That is, we treat the derivative hz = v as a new independent 
variable, so that 



(17) 



According to this formalism the Lagrangian (jl5p is replaced by the equivalent one 



L 



Auzvqq 



2nzV i) 

(1 - nir 



2nzVq 



2 + 2^,2 



-q V 



1 — nt 



+ {l-nl)q' 



— 5n1 — guz — ip{hz — v) 



where the Lagrange multiplier p ensures the equality of hz and v. The canonical momenta are p 
and Pq 

Then, partition function (IT4l) takes the form written in terms of three scalar fields q(0, w(f) and tIz{0- 



, Pv 



:dL 



z = 



J D [nz]D[v]D[q] Clip i^- L{nz,v,q)d^^ 



(18) 



OV 



(19) 



The last term in Ea.(|18p is a specific property of the Ostrogrdski methodology and we will pay a special attention 
to it because it makes the following quantum Hamiltonian a non-Hermitian one. 

Now we construct the Hamilton function H — iphz + ipvV -\- ipqq + L, which after replacing momenta by the 



corresponding differential operators: p = —i-^^, py 



and Pq = 



results in the quantum Hamiltonian: 
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It is convenient to change variables v, q, Uz to new variables r, ip, 9 connected by the relations 

T sm (p 

V = rcos(j9sin0, q — — , = cos6' (21) 

sm^ 

Then we obtain the Schrodinger equation for the quantum particle in the form 



^.2 

where 



Ho^. - ^rH. -9 casein -Scos^ 6^^ = s„*„ (22) 
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^0 = -2 + + + r + - --^^ (23) 

describes the model at the transition point at /i = 0. 

The last two terms in Ea. (P51) makes the Hamiltonian to be non-Hermitian one. Therefore, we have to consider the 
transposed counterpart of Eq. ip^ : 

H^'^n - ^rHn - g cos 0$„ - S cos^ 0$„ = e„$„ (24) 

The transposed differential operator Hq has the same form as Hq, but the sign of the last two terms in Eg . ([25)) 
is changed. This change of sign is equivalent to the change 6 — > —9, which implies that $„(r, 0, (p) — 'i'n{r,—9,(p). 
Then, the normalization condition for functions ^„ takes the form: 

— j rdr j sin ed9 J d(p'ifri{r,9,ip)-^^{r,-9,ip) = 6nm (25) 

As a result of the above manipulations the partition function Z can be considered as the partition function of 
quantum model (|22p at a 'temperature' 1/A 

Z = J2e-'^- (26) 

In the thermodynamic limit A — > cx) (A = s^^^'^T^/^L) only the lowest eigenvalue of Eq. ([22|) gives the contribution 
to Z. Thus, the free energy of the classical spin model is determined by the ground state energy Eq of the Schrodinger 
equation (12^ . The dependence of the lowest eigenvalue Eq of Eg. ip^ on the scaling parameters g, 7 and S determines 
the magnetic properties of the system. In particular, the normalized magnetization is given by 

m = — — 27 

og 



Thus, Eg. (122)) is the main result of this paper. In general, Eg. ([22| does not admit analytical solution and should 
be solved numerically. However, the limiting cases of high and low magnetic fields can be studied analytically. In the 
following we present both numerical solutions and analytical expressions for asymptotics. 

III. MAGNETIZATION CURVE AT THE TRANSITION POINT 

At first, let us consider the isotropic F-AF model at the transition point when (5 = and 7 = 0. For low magnetic 
field {g <^ 1) the ground state energy Eq can be found using the FT in g. The numerical solution of Eq.ip^ for g = 
shows that the ground state wave function does not depend on ip and 9, i.e. it satisfies the equation 

or'' r or 8 

and eo = 1.861. 

The eigenfunctions of Eg. (j22p giving the contribution to the second order in g have the form 

*n = fin (r) cos ^ sin 6* + /2„ (r) cos 9 (29) 
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where the functions /i„(r) and f2n{r) satisfy the following system of equations: 



Qj,2 J. Qj, j,2 g 



d'f2n 2df2n , 1 



■^r f2n + rfin = e„/2n (30) 



cJr^ r dr 

Normalization condition (|25p transforms for the functions /i,i(r) and f2n{i') to equation 



rdr{f2nf2m - flnflm) = ^nm (31) 



Further, we calculate the second-order correction to the ground state energy in g: 

«2 ]\.r2 

where Mon is the following matrix element 

Mon = / ■^o{r)f2n{r)rdr (33) 
Jo 

The numerical solution of Eg . ((30)) and the calculation of the sum in Ea. (l32l) gives 



Then, the magnetization m at 5 — > is 
and the zero-field susceptibility is 



£ = eo - 0.5Mg^ (34) 
■n^img + 0{g^) (35) 
1.07s8/3 



T4/3 



(36) 



Expression (1361) naturally reproduces the result found in Ref. [22| obtained by another method and confirmed by 
Monte-Carlo simulations [2^- As follows from Eg. (155)) the critical exponent of x is changed from 2 to | when a — | 
from the ferromagnetic side. 

If we assume that the hypothesis of the universality is valid, then the susceptibility x the transition point for 



the s = i F-AF chain at T is 

0-1681 

Unfortunately, the exact low-temperature asymptotic of x for s — F-AF model at a = -j is unknown. 
However, we can compare Eg. pzp with the susceptibility obtained for this model by the approximate modified spin- 
wave method (MSWT) proposed by Takahashi [3fl]. The MSWT gives x = 0.099T-'*/3 ^5, 29]. The comparison of 
MSWT result with Eg. (157)) shows that the critical exponents of both expressions are the same though the prefactors 
are different. In Ref.|29|] the transfer-matrix renormalization group (TMRG) algorithm was used for the calculation 
of the low-temperature asymptotic of x- The obtained numerical results are not fully consistent with Eq. p7|) and 
show that the exponent might actually be smaller than 4/3. However, as pointed in Ref.|29] the possible reason of 
the deviation of the TMRG results from Eg. (|37)) is that the obtainable temperatures in the TMRG calculations are 
just not low enough to observe the T~^/'^ power law predicted by Eg. (1571) . 

Now we consider the limit of large g when the magnetization is close to saturation. In this limit we expand cos 9 
near 6 = and scale the variables r and 9 as: 

Keeping in Eg. (|2^ the terms proportional to g^/** we arrive at the Schrodinger equation in a form 

2 /a^^* ^ 1 (9* ^ 1 zsinipd^^^ d'i'^x'^^ e + g^ 

\ dz^ z dz 9(^2 J J. Qj. 2 
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Fortunately, the ground state wave function and the ground state energy of Eq. (l39|) can be found exactly: 

*o(x,.,^) = Cexp(4 + f^-:^) (40) 

£0 = -g + 2.gi/4 (41) 

where C is unessential normalization constant. 

One can also calculate the next-order correction to the ground state energy (|4Tt . For this aim we estimate the 
effect of the next-order term which was omitted in Eq. (|5^ and has the form: 

9~'^ + -T- sin V'TT- (42) 



24 3 ^dip 

The calculation of the first order in perturbation gives for £0(5) the correction proportional to g^^^^: 

e, = -g + 2g'/' + ^ (43) 
Then, the asymptotics for the magnetization and the susceptibility for h ^ T"^/^ are 



3T 



8/l7/4sl/4 



2/j3/4g5/4 V/l3/2s5/2 



-12 



(45) 



It is interesting to compare the leading terms of this expansion with the spin- wave expansion of the magnetization for 
the spin-s quantum F-AF chain at a = ^. We have checked that this expansion reproduces Eq. pi)) . The second term 
in Eq. (j44p corresponds to the result of the linear spin- wave theory, but the third one includes spin- wave interaction 
effect and, therefore, the coincidence is not trivial. Certainly, we can not prove that both expansions coincide in all 
orders in small parameter g~^^'^. Nevertheless, the coincidence of the leading terms of m{g) for the quantum and the 
classical model gives a promise that the universality is valid at the transition point of the F-AF model. 

We complete this subsection with the results for the spin correlation function. It can be shown [l^] that the spin 
correlation function {SziO)Sz{l)) has the form 

{SMSzil)) = J2 (*o|*n)' exp[-ri/3(e„ - £0)^] (46) 

As follows from Eq.(j46| the spin correlation function exponentially decays on long distances I ^ T^^^^, and the 
correlation length is governed by the lowest eigenstates of Eg. ([2^ . In the case of absence of the magnetic field {g = 0) 
all the eigenvalues £„ are real and several lowest levels was calculated in [22], which gives the asymptotic for the 
correlation length = 1.04r^^/^ at 5 0. 

In the high field limit ((/ ^ 1) there are three lowest excited states having equal real part of their eigenvalues: 

ei = -.9 + 4.gi/4 + 0(5-1/2) 
£2,3 = -.9 + 2gi/4(2±,) + 0(^-1/2) (47) 

According to Eq. pS)) the presence of the imaginary part in eigenvalues (|T7)) causes the oscillations on the background 
of the exponential decay of the correlation function. The imaginary part of the eigenvalues determines the period of 
the oscillations while the real part determines the correlation length. According to Eas. (|15)) and (|47p the asymptotic 
of the correlation length in the high field limit is 

^0 = 2.5/12/,l/4 (48) 

So we see that the correlation length is defined by the temperature for h <^ T^/^ and by the magnetic field when 
h 3> T**/^. We note that the ratio of the correlation lengths in these limits is proportional to g^^^- The crossover 
between these two regimes occurs at g ~ 1. 

In general, the solution of Eq.(|22|) and the computation of £o(.g) and m{g) has been obtained numerically. The 
dependence m{g) at the transition point is shown by thick solid line in FigHJ 
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FIG. 1: Magnetization curves m{g) at 7 = for some values of scaled anisotropy 5. The case 5 = corresponds to the transition 
point. 



IV. MAGNETIZATION CURVE IN THE HELICAL PHASE 



In this section we consider the behavior of the magnetization in the helical phase in the vicinity of the transition 
point, when a > j, but the anisotropy is zero 3 = 0. For a > j the ground state has the helical type of long-range 
order (LRO) with the wave-number kh = cos^^(1/4q;). The saturation field hs at a close to j is hs = s(4q! — 1)^. 
The ground state magnetization is given by to = h/hg for h < hg and to = 1 for h > h^. At finite temperature the 
helical LRO is destroyed by thermal fluctuations and thermodynamic quantities have singular behavior at T — >■ 0. 

The behavior of the system in the case of absence of the magnetic field was studied in detail in Ref. [1^ . It was 
shown that with the increase of the temperature the gapless excitations over the helical ground states (spin waves) 
smear the (5-peaks of the static structure factor at A: = ±kh and shift the peaks to A: = 0. Finally, at Tc = 0.9257'^/^ 
the maximum of the spin structure factor reaches fc = defining the Lifschitz boundary, so that the helical type of 
the spin correlations for T > Tc disappears. 

Most likely, the hypothesis of the universality of the function m{g) breaks down for a > |; because the excitations 
above the ground state are different in the quantum and in the classical F-AF chain. Nevertheless, as was shown in 
Ref. [i?! some peculiarities of the low-temperature behavior of the classical model at a > -j is qualitatively similar to 
that for the quantum s = ^ chain. For example, the temperature dependence of the zero-field susceptibility is in a 
qualitative agreement with the numerical data for the quantum s = ^ model and is in accord with the experimental 
data for the real edge-sharing compounds. 

The finding of the magnetization in the helical phase reduces to the solution of Eq.([22|) for 7 > and can be 
analyzed in full analogy with the case 7 = 0. For low magnetic field {g — )■ 0) the magnetization is to ^ g and can be 
represented as 

m = Ag(7) (49) 
hs 

The function G(7) is found from the solution of Eg. ([50]) where the terms —\"yr'^fi(2)n are added to the first (the 
second) equation of (|30p . In fact, G(7) coincides with the normalized zero-field susceptibility obtained before in 
Ref.[24|. Therefore, we do not present this function here. We note only that G(7) vanishes at 7 — > 00, tends to finite 
value at 7 — >■ and has a maximum at 7 ~ 2.2. 

In high magnetic field limit we use rescaling (P5| for Eq.([22) and keep the leading terms. The obtained equation 
repeats Eg. ([5^ with the additional term —^^/hjhz^'^. The ground state wave function of this equation has the 
form similar to Eg. pO)) : 

* = Gcxp ^-az^ + ^2;a;cos(/3 - 2ax^^ (50) 

and the ground state energy is 



£0 = -g + 8agi/^ 



(51) 




FIG. 2: Magnetization as a function of h/hs for several values of parameter 7 for isotropic case [5 = 0). 



where 

a= i^l- (52) 
Eg . ([50)) is valid for high fields when h> hg. The asymptotic of the magnetization curve in this limit has the form 

m = 1 ^ ' ' (53) 

As follows from Eg . ((531) the temperature-dependent correction to to = 1 at /i ^ ft.^ is proportional to 7^"^/^. 

The magnetization curves for several values of 7 as a function of hjhs are shown in Fig 12] together with the ground 
state magnetization (7 — >■ cxi). The zero-field susceptibility defines the slopes of the magnetization curves for small 
hjhs and as follows from Figl5]this slope can be both larger and smaller than the ground state value. Such behavior 
of the magnetization follows from the non-monotonic dependence of G on 7. 



V. EASY-AXIS ANISOTROPY 



Up to now we considered the isotropic F-AF chain. At the same time it is important to study the influence of the 
anisotropy of exchange interactions on the low-temperature thermodynamics. In this subsection we pay our attention 
mainly to the dependence of the zero-field susceptibility on the anisotropy at a = i. 

At first, we briefly review the effect of the anisotropy on the susceptibility for the classical and quantum ferromag- 
netic model (a = 0). At T — >■ and in the scaling limit the magnetization of the classical ferromagnetic chain is given 
by Eg. ip?)) . where Eq is the lowest eigenvalue of the equation pisl [T9|: 

1 \ 

- 7: ~ 7: cot 6'— - (5f cos^ 0^ - gp cos 6** = e'^ (54) 

2 dO'' 2 dO 

In this equation Sp = {A ~ l)s^/T'^ and gp — s^h/T'^ are the scaling parameters. 

To find the susceptibility in the limit 5^ <C 1 and gp <^ 1 we can use the PT mV = —Spcos^O — gpcosO. At 
y = the eigenfunctions of Ea. (l54|) are the Legendre polynomials Pi{cos9) with the eigenvalues £1 — l{l + l)/2. The 
PT in the lowest orders in V gives 

_ Sp gl AAg^Sp 

Then, according to Ea. (l?fl) the zero-field susceptibility in the limit of weak anisotropy 5f ^ (A — 1 ^ T^) is 

2s4 / 44s4(A~l) 



1 + ^1^) (56) 



In the opposed limit 5p ^ 00 and gp = Eq. (|5^ has two almost degenerated lowest eigenvalues corresponding to 
the states with even and odd parity with respect to exchange 6 ^ {tt—9). The tunnel splitting between these states can 
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be found with the exponential accuracy in the WKB approximation [3i|: Ae ~ exp(— 2a/27p). The term {—gp cosO) 
in Eq.(|54|) has non-zero matrix element between the states with even and odd parities, so that the contribution to 
the second order PT in gj? is given by 

eo- -5|exp(2V2^) (57) 

and the zero-field susceptibility is 

X-^exp(^) (58) 



where AE = 2s V2(A - 1). 

As follows from Ea. (l58l) the susceptibility diverges exponentially at T — >■ and the value of the thermal gap AE 
is equal to the kink energy (or one-half of the energy of large soliton) of the weakly-anisotropic ferromagnetic chain. 
As it is known (s^ the soliton energy coincides with the energy of the multimagnon bound states of the quantum 
anisotropic ferromagnetic chain. It is interesting to note that the susceptibility of the easy-axis anisotropic s ~ ^ 
ferromagnetic chain found on the base of the Gaudin formalism in Ref.[33] behaves at T — >■ as x ~ exp{AE/T)/T 
and AE is the same as given in Eg. ([55)) . This fact manifests the close relation between the magnetic properties of the 
quantum and the classical anisotropic ferromagnetic chains. 

We will show that this resemblance remains in the F-AF model at the critical point a = j. For example, it was 
shown by us in Refs.[25l |35| that the energies of multimagnon bound states in the quantum model and the energy of 
large classical solitons at a = j are both proportional to (A — 1)'^/'* though the numerical coefficients are different. We 
have also shown [2^ that the susceptibility of the quantum s = i F-AF model diverges at T — ?> exponentially, i.e. 
X ~ exp(Ai?/T) and the thermal gap AE equals one-half of the energy of the multimagnon complexes. We will show 
below that the susceptibility of the classical model at (5 ^ 1 has similar exponential behavior and the corresponding 
thermal gap is the classical kink energy. 

Similar to the pure ferromagnetic chain the calculation of the zero-field susceptibility reduces to the computation 
of the ground state energy of Ea. (p2)l in the second order mV = —gcosd. Then the susceptibility can be represented 
as X = 2(s7r)4/3/(^) where 

m = (59) 

and 'i/n and e„ are the eigenfunctions and the eigenvalues of Ea. (l22t at g = 0. 

It is convenient to introduce the normalized susceptibility x — (A — l)x and the normalized temperature T = 
T/s^(A — 1)^/"' (T — (5^'^/^), so that x is a function of T only. The function x{T) can be found explicitly in the 
limits of large and small values of T (at small and large 6 correspondingly) . For high temperatures T ^ 1 the states 
giving the contributions to the sum in Eq. ((59|) are separated from the ground state by finite gap and the numerical 
calculation of this sum gives 

X = 1.077'-'*/^ -I- 0.87T'-^/3 (60) 

The calculation x(r) for small T is more complicated. At T — > ((5 — > oo) we can expand the term (5cos^ 9 up to 
9'^ in the Schrodinger equation (l22t . Then the ground state wave function has the form similar to Eq. pO)) and the 



energy Eq = —T -I- I^I^T . Further, we can compute perturbative corrections to Eq from omitted anharmonic 
terms to obtain the ground state energy in a form 

i?o=eo + ^a„f"+2/3 (61) 



If we expand the term 8 cos^ 9 near the second minimum 9 = tt we would obtain the result identical to Eq. (1611) , so 
that we have two degenerated states. But there is a non-perturbative tunnel splitting which is not captured by the PT. 
The splitting is exponentially small at T as demonstrated in Figl3] On the other hand, these quasi-degenerated 
states have a non-zero matrix element in Eg. (1591) and, therefore, the splitting between them determines the behavior 
of X at f -> 0. 

The most convenient way to evaluate this splitting is the calculation of the original functional integral (fH)) . Cer- 
tainly, the exact calculation of this integral is impossible and, therefore, we use a semiclassical approximation. In this 
approximation the functional integral (|14p is represented as the sum of the contributions of the classical paths in Z 



11 








— even 

— odd 



1.5 



2.5 



FIG. 3: Dependence of two lowest eigenvalues of Eq. (|22P on the normalized temperature T — T/s^(A — I)"'/'* at 7 = g = 0. 






FIG. 4: Instanton solution of Ea. (|62[l with boundary conditions: 0{ — oo) — 0, 6l(oo) — n. 



minimizing the Euclidean action and the paths which are close to the classical ones. The minimization of the energy 
fmictional W{6, ip) (Eq. ([T5|) ) gives = and the following Euler equation for 6{£): 



1 



r4/3 



sin(26l) = 



(62) 



We note that there are a few classical solutions of this Euler equation. Two of them 6 = {) and 6 — tt correspond 
to trivial ferromagnetic configurations. A systematic expansion around these saddle points is equivalent to a purely 
perturbation expansion of the eigenvalues of the Schrodinger equation (1221) . The corresponding contribution to the 
partition function is Z'"' = exp(— AE'o) where Eq is given by Eg. dBTj) and A is the scaled system length pT|) . The 
tunnel splitting is given by an instanton contribution to the functional integral ([H]) (33 | . The classical solution of the 
Euler equation corresponding to the instanton satisfies the boundary condition: 9(—oo) = 0, 0(oo) = tt or 6(—oo) = tt, 
9(00) = 0. This solution of the Euler equation can be found numerically. Actually, it coincides with the solution for 
the kink excitation in the F-AF chain found by us in Ref.js^. The dependence in this solution is shown in FigH) 
According to the result of Ref.[35|, the classical action corresponding to instanton is Sq = 2/T. The summation of 
the contributions to Z using the semiclassical approximation can be performed by a standard way [33 |. As a result, 
the partition function is represented in a form 



(63) 



where Z^"^ is defined above, Z^^^ is the instanton-antiinstanton contribution (lA), Z'^^^ is the lAIA contribution and 
so on. 

Using the usual approximation of the semiclassical method (in particular, neglecting instanton-instanton interac- 
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FIG. 5: Dependence of the normalized magnetic susceptibility x = (A — \)x on the normalized temperature T — T/s^(A — 1)'''''* 
(thick solid line). Dashed and thin solid lines are asymptotics of x at T ^> oo and T — > given by Eqs. (|60|I and 
correspondingly. 



tions) we have 



^(2) = 



-XEq 



-\Eo 



and so on (here A is the scaled system length ([TT]) ). 
Summing (|63| we arrive at 

Z = e"^-^" cosh(Ae"-^°) 
On the other hand, the partition function at A — >■ cx) can be represented as 



(64) 



(65) 



(66) 



where Es and Ea are the energies of the lowest states with even and odd parity with respect to exchange (r, 6) ^ 
(-r, TT - 9). Then, Es ^ Eq- e'^" and -Ba = Eq + e~^°. The tunnel splitting is 



AE = Es- Ea = 2exp(-2/f') 
Then the susceptibility at T — >■ to the exponential accuracy is given by: 

"2s2(A- 1)3/4" 



X ^ exp 



T 



(67) 



(68) 



The thermal gap in Eg. ([68)) is the kink energy of the weakly anisotropic classical F-AF chain. It is interesting to 
compare x fEq. ((65|) ') with the susceptibility of the quantum F-AF model at a = | ^sj. The susceptibility for both 
models shows the exponential dependence with the thermal gap proportional to (A — l)^/^. If we use Eg. (155)1 for 
s = 2 case we find that the thermal gap is i(A - 1)3/4 ^Yiile in fact it is 0.35(A - i f/^ [H, i.e. the numerical 
coefficients at (A — 1)^/4 are slightly different. 

Eqs.(|Bn)) and (1551) give asymptotics of x{T) for small and large values of T. In general case the function xi^T) has 
been calculated numerically and the dependence xC^) is shown in FiglS] together with asymptotics of x for small 
and large T. As it can be seen from Figl5] the dependence x on T is characterized by two types of behavior: x is 
proportional to T'~4/3 at T 3> 1 and grows exponentially at T — > 0. The crossover between two regimes occurs at 
f - 1. 

The calculation of the susceptibility of the anisotropic F-AF model can be expanded to the case 7 > 0. We do 
not dwell on details of these calculations. We notice only that x becomes the function of T and of a parameter 
fi — (4a — A — 1. x(/i,T) as a function of T has a minimum for /i > /xq ~ 1 and is finite at T — >■ 0. The 
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normalized susceptibility x diverges at T — >■ for /i < /iq. The line [Aa — 1) = /iq %/ A — 1 can be identified with the 
boundary between the ferromagnetic and the helical phase. 

Finally we give the results for the magnetization of the anisotropic F-AF model at a = i. According to Eg s . (160]) 
and (|68p m{g) at g — > is 

m{g) = 1.07g(l + 0.81(5), (5 < 1 

m{g) - 5exp(2(5^/''), (5 > 1 (69) 

The behavior of the magnetization in the high magnetic field limit is obtained by analogy with that for the isotropic 
case (see Eq. (|44| \ Then m{g) at ,g > 1 is 

"^-^-2(^T^ + ^((^ + 

The magnetization curves m{g) for several values of 5 obtained by the numerical solution of Eg. ip^ is shown in 
FigUl Its behavior in the limits g <C 1 and 5^1 agrees with Eas. ((M)) and ((7U)) . 



VI. CONCLUSIONS 



We have studied the low-temperature magnetic properties of the classical anisotropic F-AF chain in the vicinity of 
the transition point from the ferromagnetic to the helical ground state. This means that the frustration parameter 
a = J2/I Ji| is close to its critical value a = \ and the anisotropy of the exchange interaction (A — 1) is weak. In the 
vicinity of the transition point the nearest spins in the ground state are directed almost (or even exactly) parallel to 
each other. Therefore, in the low-temperature limit when the thermal fluctuations are weak, we can use the continuum 
approximation and represent the partition function as a functional integral over the spin vector field. In the obtained 
energy functional the model parameters (a — |;), (A — 1) and the magnetic field h are scaled by the temperature and 
form three independent scaling parameters 7, 5 and g defined in Eq. (jl2p . This implies that we considered the scaling 
limit when a— J-^, A— >-l, ft,— >-0 and T — >■ 0, but the values of the scaling parameters 7, 5 and g are finite and 
govern the low-temperature thermodynamics of the F-AF model near the transition point. 

The derived functional integral for the partition function was treated as a path integral of the quantum mechanics. 
The peculiarity of this path integral is that the Lagrangian contains the second order derivative. To handle with 
this problem we used the special Ostrogradski prescription, which allowed us to obtain the quantum Hamiltonian 
corresponding to such path integral in a special unusual form. Then the dependence of the lowest eigenvalue of the 
Hamiltonian on the scaling parameters determines the magnetization curves of the system. The eigenvalue problem 
has been solved numerically and explicit expressions for the magnetization was obtained in the limits of low and high 
magnetic fields. 

It is known [18] that the magnetization curve for the pure ferromagnetic chain has a universal form when plotted 
against the scaled magnetic field gp = s^h/T'^, and this curve is valid for any value of spin s including the classical 
limit s — > cx). We suppose that such universality remains for F-AF model at the transition point against the scaling 
parameter g = hs^^'^ /T^^^ . If this is the case the obtained magnetization curves for the classical model can be easily 
recalculated to the quantum spin case. To validate this hypothesis one needs to compare the obtained classical results 
with the magnetization of the s = ^ F-AF model. Unfortunately, the exact thermodynamics of the latter model is 
unknown. Nevertheless, there are two indirect arguments supporting this conjecture. First is that the obtained critical 
exponent | in the temperature dependence of susceptibility at the transition point coincides with that obtained in the 
MSWT method. The second argument is that three leading terms of the spin-wave expansion of the magnetization of 
the quantum model coincide with those for the classical model. Certainly these two facts do not prove the proposed 
hypothesis and the question about its validity remains open poj . In this respect the numerical calculations of the 
magnetization as a function of the magnetic field at T — ^ for s = ^ and s = 1 are very desirable. 

Probably, the hypothesis of the universality of the function m{g) (if any) breaks down for a > ^ because the 
excitations above the ground state are different in the quantum and in the classical F-AF chain. Nevertheless, as was 
shown in Ref.l22| some peculiarities of the low-temperature behavior of the classical model at a > |; is qualitatively 
similar to that for the quantum s — ^ chain. For example, the temperature dependence of the zero-field susceptibility 
is in a qualitative agreement with the numerical data for the quantum s = ^ model and is in accord with the 
experimental data for the real edge-sharing compounds. 

We have studied the influence of the easy-axis anisotropy on the behavior of the susceptibility at the transition 
point. It is shown that even weak anisotropy essentially changes x- In the low-temperature limit the susceptibility 
diverges exponentially in contrast with the isotropic case where the divergence is of a power-like type. We note 
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that such behavior of the susceptibihty takes place in the quantum s = | F-AF chain f25| and the corresponding 
thermal gap has the same functional form as the classical one. This fact confirms the close relation between the 
low-temperature magnetic properties of the quantum and classical F-AF model in the ferromagnetic part of the phase 
diagram. 
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